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Abstract 

This paper is concerned with giving the proof that there is a gen- 
eral decouphng property of vacuum and nonvacuum Einstein equa- 
tions written in variables adapted to nonholonomic 2-1-2 splitting. We 
show how such a geometric techniques can be applied for construct- 
ing generic off-diagonal exact solutions of Einstein- Yang-Mills-Higgs 
(EYMH) equations. The corresponding classes of solutions are de- 
termined by generating and integration functions which depend, in 
general, on all space and time coordinates and may possess, or not. 
Killing symmetries. The initial data sets for the Cauchy problem and 
their global properties are analyzed. There are formulated the criteria 
of evolution with spacetime splitting and decoupling of fundamental 
field equations. Examples of exact solutions defining black ellipsoid 
and solitonic configurations are provided. 
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1 Introduction 



The issue of constructing exact and approximate solutions of gravita- 
tional and matter field equations is of interest in mathematical relativity, 
particle physics and for various applications in modern cosmology and as- 
trophysics. The Einstein, Yang-Mills, Higgs and other fundamental field 
equations in physics are sophisticate systems of nonlinear partial differen- 
tial equations (PDE) which are very difficult to be integrated and studied 
in general forms. There were elaborated various geometric, analytic and 
numeric methods, see PmSlEHHElElEllS] and references therein. 

In this work, we address again and develop a geometric approach (the 
so-called anholonomic deformation method, see [9l [T0\ lllj ) to construct- 
ing exact solutions of PDE describing generic off-diagonal and nonlinear 
gravitational, gauge and scalar field interactions. The goal is to provide 
new results on decoupling propert}!!] of Einstein- Yang-Mills-Higgs, EYMH, 
equations and integration of such systems. Certain issues on the Cauchy 
problem and decoupling of Einstein equations and the initial data sets and 
nonholonomic evolution (with non-integrable constraints and/or with re- 
spect to anholonomic frames) will be analyzed. Examples of off-diagonal 
solutions for EYMH black holes/ellpsoids and solitonic configurations will 
be given. 

In Section [2] we provide the main Theorems on decoupling the Einstein 
equations for generic off-diagonal metrics with one Killing symmetry and 
consider extensions to classes of "non-Killing" solutions with coefficients 
depending on all set of four coordinates. We prove also that the decoupling 
property also holds true for certain classes of nonholonomic EYMH systems. 
The theorems on generating off-diagonal solutions are considered in Section 
[3l Section H] contains a study of the Cauchy problem in connection to the 
decoupling property of Einstein equations. Next two sections are devoted to 
examples of generic off-diagonal exact solutions of the EYMH equations. In 
Section [5] there are constructed nonholonomic YMH vacuum deformations 
of black holes. Ellipsoid-solitonic non-Abelian configurations are analyzed 
in Section [6l Appendix |X] contains a survey on the geometry of nonholo- 
nomic 2+2 splitting of Lorentz manifolds. The most important formulas 
and computations which are necessary to prove the decoupling property of 
Einstein equations are given in Appendix [Bj 

^it is used also the term "separation" of equations, which should not be confused with 
separation of variables 
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2 Decoupling Property of Einstein Equations 



In general relativity (GR), the curved spacetime (Kg) is defined by a 
pseudo-Riemannian manifold V endowed with a Lorentzian metric g as a 
solution of Einstein equations H The space of classical physical events is mod- 
elled as a Lorenzian four dimensional, 4-d, manifold V (of necessary smooth 
class, Housdorf and paracompact one) when the symmetric 2~covariant ten- 
sor g = {gap} defines in each point u € V and nondegenerate bilinear form 
on the tangent space T^V, for instance, of signature (+ -|- H — ). The as- 
sumption that TuV has its prototype (local fiber) the Minkowski space M'^'^ 
leads at a causal character of positive/negative/null vectors on V, i.e. for 
the module X of vectors fields X,Y, ... E ^iV), which is similar to that in 
special relativity. 

Let us denote by ea and a local frame and, respectively, its dual frame 
[we can consider orthonormal (co) bases], where Greek indices a,/3, .. may 
be abstract ones, or running values 1, 2, 3, 4. For a coordinate base u = {u°'} 
on a chart U C V, we can write Ca = da = d/du^ and = du^ and, for 
instance, define the coefficients of a vector X and a metric g, respectively, 
in the forms X = X^Ca and 

g = 5a/3(w)e"0e'^, (1) 

where gaj3 ■= g(ea,e^)|l We consider bases with non-integrable (equiva- 
lently, nonholonomic/anholonomic) 2 + 2 splitting for conventional, hori- 
zontal, h, and vertical, v, decomposition, when for the tangent bundle TV 
:= 1J„ TuV a Whitney sum 

N : TV = hV®vV (2) 

is globally defined. Such a nonholonomic distribution is determined locally 
by its coefficients N^(u), when N =N^{x, y)dx^®d/dy°', where ti" = (x*, y°) 
splits into h-coordinates, x = (x*), and v-coordinates, y = (y°), with indices 
running respectively values i, k, ... = 1,2 and a, 6, c, ... = 3, 4o 

^We assume that readers are familiar with basic concepts and resuhs on mathematical 
relativity and methods of constructing exact solutions outlined, for instance, in above 
mentioned monographs and reviews. 

^The summation rule on repeating low-up indices will be applied if the contrary will 
be not stated. 

*We note that the 2+2 splitting can be considered as an alternative to the well known 
3+1 splitting. The first one is convenient, for instance, for constructing generic off- 
diagonal solutions and elaborating models of deformation and/or A-brane quantization 
of gravity, but the second one is more important for canonical/loop quantization etc. 
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Any spacetime {V, g) can be equipped with a non-integrable fibred struc- 
ture ([2]) and such a manifold is cahed nonholonomic (equivalently, N-anholo- 
nomic). We use "boldface" letters in order to emphasize that certain spaces 
and geometric objects/constructions are "N-adapted", i. e. adapted to 
a h-v-splitting. The geometric objects are called distinguished (in brief, 
d-objects, d-vectors, d-tensors etc). For instance, we write a d-vector as 
X = {hX,vX) for a nonholonomic Lorentz manifold/spacetime (V,g). 

On a spacetime (V,g) , we can perform/adapt the geometric construc- 
tions using " N-elongated" local bases ( partial derivatives), e^, = (ej,ea), 
and cobases (differentials), = (e*,e"), when 

e, = d/dx' - N^{u)d/dy^ ea = da = d/dy\ (3) 
and = dx\ e" = dy" + N,^{u)dx\ (4) 

Such (co) frame structures depend linearly on N-connection coefficients be- 
ing, in general, nonholonomic. For instance, the basic vectors ([3]) satisfy 
certain nontrivial nonholonomy relations 

[e„, e/j] = BaB/s - e^e„ = Wj^e^, (5) 

with (antisymmetric) nontrivial anholonomy coefficients 

= daNl Wf, = n% = ej {Nn - e.(7V;). (6) 

Any spacetime metric g = {gajs} via frame/coordinate transforms 
can be represented equivalently in N-adapted form as a d-metric 

g= 9ij{x,y) (^e^ + gab{x,y) e"- ®e^, (7) 

or, with respect to a coordinate local cobasis du" = {dx'^,dy°'), as an off- 
diagonal metric 



where 

gij + N^N!lgab N^gae 
Nfgte gab 



(8) 



A metric g is generic off-diagonal if ([8]) can not be diagonalized via coordi- 
nate transforms. Ansatz of this type are used in Kaluza-Klein gravity when 
N^{x,y) = V'^{x)y"' and are " compactified" extra-dimensions coordi- 
nates, or in Finsler gravity theories, see details in \12\ 111]. In this work, we 
restrict our considerations only to the 4-d Einstein gravity theory. The prin- 
ciple of general covariance in GR, allows us to consider any frame/coordinate 



4 



transforms and write a spacetime metric g equivalently in any above form 
([1]), (l8|) and/or d?]). The last mentioned parametrization will allow us to 
prove a very important property of decoupling of Einstein equations with 
respect to N-adapted bases ([3]) and (jl]). 

Via frame/ coordinate transforms Bq. = e°'a{x,y)ea' , g^^ = e^^e^ ^a'/j" 
a metric g d?]) can be written in a form with separation of v-coordinates 
and nontrivial vertical conformal transforms, 

g = gidx'' (g) dx^ + uj'^hahaG'^ ® e", (9) 
= dy^ + {wi + w_j) dx^, e'^ = dy'^ + (rij + nj dx', 

were 

gi = gi{x^),ga = ^'^{^\yl ha{x'',y^)h^{x^,y^), 

Nf = Wi{x\y'')+Wi{x\y^),Nf = ni{x\y^)+ni{x\y^), (10) 

are functions of necessary smooth class which will be defined in a form to 
generate solutions of Einstein equations!! 

The aim of this section is to prove that the gravitational field equa- 
tions in GR, in the vacuum cases and certain vary general classes of matter 
field sources decouple for parametrizations of metrics in the form (jlOp . For 
convenience, we present in Appendix |A] a brief review on the geometry of 
nonholonomic 2+2 splitting of Lorentz manifolds and Einstein equations. 

2.1 Splitting of (non) vacuum gravitational field equations 

We shall use brief denotations for partial derivatives, a* = da/dx^,a' = 
da/dx'^,a* = da/dy^,a° = da/dy'^. The equations will be written with 
respect to N-adapted frames of type ([3]) and (jH). 

2.1.1 Off diagonal spacetimes with Killing symmetry 

We use an ansatz ([9]) when f^ = l,/i3 = l,wij = and = in data 
([To]) and T = for (jA.lSp . Such a generic off-diagonal metric does not 
depend on variable y^, i.e. d/dy^ is a Killing vector, if /14 = 1. Nevertheless, 
the decoupling property can proven for the same assumptions but arbitrary 
hji{x^,y'^) with nontrivial dependence on y^. We call this class of metrics 
to be with weak Killing symmetry because they result in systems of PDE 
()A.12p as for the Killing case but there are differences in ()A.13P if /14 7^ 1. 

^ There is not summation on repeating "low" indices a in formulas (|10[1 but such a 
summation is considered for crossing " up-low" indices i and a in ([U ) . We shall underline 
a function if it positively depends on but not on tf^ and write, for instance, n^ix'' , y'^). 
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Theorem 2.1 The Einstein eqs IIA.1!3\) and iA.13\) for a metric g U0\) with 

= = 1 and Wi = n^=0 and T = in matter source T^^ IIA.15\) are 
equivalent, respectively, to 



Rl = i?^ = ;T^[5r-#-¥^+ffi'-#-¥^] = -"T, (11) 



and 

w* = ((9i - u;i)ln|/i4|,5feWi = (9iWfc, (15) 
nklu = t^feZLi, fi* = 0, dtUk = dkn^. 

Proof. See Appendix [Bj □ 

Let us discuss the decouphng (sphtting) property of the Einstein equa- 
tions with respect to certain classes of N-adapted frames which is contained 
in the system of PDE pip ~ (|15p . For instance, the first equation is for a 
2-d metric which always can be diagonalized, [51,52], and/or made to be 
conformally flat. Prescribing a function gi and source "T, we can find 52 > 
or inversely. The equation (|12p contains only the first and second deriva- 
tives on d/dy^ and relates two functions h-^ and /i4. Prescribing one of such 
functions and source T, we can define the second one taking, respectively, 
one or two derivations on y^. The equation (I13p is a linear algebraic system 
for Wk if the coefficients ha have been already defined as a solution of p2p . 
Nevertheless, we have to solve a system of first order PDE on and y^ in 
order to find resulting in zero torsion conditions (jlSp . Such conditions 
do not allow a complete decoupling because the first equations relate Wi 
to H = ln|/i4| via corresponding first order PDE. Nevertheless, it is possi- 
ble, for instance, to integrate such solutions for any prescribed H (see more 
details below, in Remark-Example 13. ip . The forth equations (|14p became 
trivial for any n* = if we wont to satisfy completely such zero torsion 
conditions If] A nontrivial function h^ is explicitly present in the conditions 
([15]). If such restrictions are satisfied, this allows us to eliminate h^ from 
(see details in the proof of above theorem). 



^Nontrivial solutions and nonzero torsion configurations present interest in modified 
theories of gravity, see such examples in Ref. [IT] 
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We conclude that the Einstein equations for metrics with one Kilhng 
symmetry can be such way parametrized with respect to N-adapted frames 
that they decouple and separate into "quite simple" PDE for h-components, 
gi, and then for v-components, ha- The N~connection coefficients also sep- 
arate and can be defined from corresponding algebraic and/or first order 
PDE. The "zero torsion" conditions impose certain additional constraints 
(as some simple first order PDE with possible separation of variables) on 
N-coefficients and coefficients of v-metric. 

In a similar form, we can decouple the Einstein equations for spacetimes 
with one Killing symmetry on d/dy^. 



Corollary 2.1 The Einstein eqs IIA.1S\) and IIA.13\) for a metric g [T0\) with 
uj = h/i = 1 and Wi = Ui = 
equivalent, respectively, to 



(J = /i4 = 1 and vui = Ui = and T = in matter source T ^ IIA.15\) are 



3l 32 If .9*.92 (.92) , // 9l92 (.9l)^i vy- 

^1 ^2 = [92 - -7, 7. — + 9i - -7 7; — \ = - JL , (16) 

2gi52 2gi 2g2 2g2 2gi 

^ Zife r 00 (Z^-s) h^lki^ ( dkh.3 dkh4\ dk^s 



and n° = {d^ - In \h^\, {dk - n^^n, = {di - n^)n^, (20) 
w^hl = dkh3,w° = 0,diW^ = dkW^. 

Proof. It is similar to that for Theorem 12.11 provided in Appendix [Bj 
We do not repeat such computations. □ 

Using above Theorem and Corollary and mutual transforms of the sys- 
tems of equations, we can formulate: 

Conclusion 2.1 The nonlinear systems of PDE corresponding to Einstein 
equations iA.l^} and iA.13\) for metrics g il(J\) with Killing symmetry on 
d/dy4, when a; = /13 = 1 and u^j = Hj = and T = m matter source T^^ 
/iA.15\) . can be transformed into respective systems of PDE for data with 
Killing symmetry on d/dy^, when w = /14 = 1 and Wi = Ui = and T = 0, 
if h3{x\y^) luix^y*)^ h4{x\y^) h^ix^y^^), Wk{x\y^) n^ix^y'^) 
and nk{x'^ ,y^) — )■ Wi^{x^,y^)- 
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The above presented method of nonholonomic deformations can be used 
for decouphng the Einstein equations even some metrics do not possess, in 
general, any Kilhng symmetries. The generic nonUnear character of such 
systems of PDE does not allow us to use a principle of superposition of 
solutions. Nevertheless, certain classes of conformal transforms for the v- 
components of d-metrics and nonholonomic constraints of integral varieties 
give us the possibility to extend the anholonomic deformation method to 
"non-Killing" vacuum and nonvacuum gravitational interactions. In next 
two subsections, we analyze two possibilities to decouple the Einstein equa- 
tions for metrics with coefficients depending on all spacetime coordinates. 

2.1.2 Preserving decoupling under v conformal transforms 

This property is stated by 

Lemma 2.1 The Einstein equations IIA.1!3\) for geometric data iB.l]} . i.e. 

the system of PDE ^2\) - [14\ ), do not change under a "vertical" conformal 
transform with nontrivial Lo{x^,y"') to a d-metric I110\) if there are satisfied 
the conditions 

d^oj — WiUJ* — niU}° = and T^^ = 0. (21) 

Proof. It follows from straightforward computations when coefficients 
9iix''),g3 = /i3(a;^2/^),54 = /i4(x^y^)/l4(x^,y^),iV^^ = Wi{x'',y^),N^ = 
ni{x^,y^)) are generalized to a nontrivial ui{x^,y"') with '^g^ = oj'^h^ and 
'^54 = Lo'^hih^. Using respectively formulas (|A.2p . (jA.3p . ()A.6p and (|A.7p . we 
get distortion relations for the Ricci d-tensors (|A.lip . 

"i^'^fe = R\ + and ^Ru = Ru = 0, 

where i?"^ and R\,i are those computed for a; = 1, i.e. p^ - p^ . The values 
and Z\ are defined by a nontrivial io and computed using the same 
formulas. We do not repeat here such details provided in Refs. [lOl [H] for 
/i4 = 1 because a nontrivial hj^ does not modify substantially the proof that 
'^Z'^j, = if the conditions are satisfied. □ 

Using Theorem 12.11 Corollary [2T] Conclusion 12. 11 Lemma |2. 11 we prove 

Theorem 2.2 A d-metric 

g = 5j(2;^)(ix*(g)dx* + a;^(x'',y'') (/i3e^(g)e^ + /i4 /i4e'^(g)e'^) , 
e^ = dy^ ^ Wiix^ ,y^)dx\ = dy^ ^ ni{x^)dx\ (22) 
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satisfying the PDE |77]j- [T5\) and dkco — wiuj* — niUj° =0, or a d-metric 

= dy^ + Wi{x'')dx\ = dy^ + ni{x'',y^)dx\ (23) 

satisfying the PDE [TS\} - (2U\) and dkOJ — w.iUJ* — n^uj° = 0, define, in general, 
two different classes of generic off-diagonal solutions of Einstein equations 
fA.l^} and fA.13\} with respective sources of type ^A.15\} . 

Both ansatz of type ()22p and (j23p consist particular cases of parametriza- 
tions of metrics in the form ([9]). Via frame/coordinate transform into a finite 
region of a point ''n G V any spacetime metric in GR can be represented 
in an above mentioned d~metric form. If only one of coordinates y"" is time- 
like, the solutions of type (j22p and (|23p can not be transformed mutually 
via nonholonomic frame deformations preserving causality. 

2.1.3 Decoupling with effective linearization of Ricci tensors 

The explicit form of field equations for vacuum and nonvacuum gravi- 
tational interactions depends on the type of frames and coordinate systems 
we consider for decoupling such PDE. We can split such systems for more 
general parameterizations (than ansatz (j22p and (I23p ) in a form ([9]) with 
nontrivial cj and N-coefficients in (jlOp . This is possible in any open region 
[/ C V where for computing the N-adapted coefficients of the Riemann and 
Ricci d-tensors, see formulas ()A.6P and ()A.7p . we can neglect contributions 
from quadratic terms of type T ■ T but preserve values of type S^F. For 
such constructions, we have to introduce a class of N-adapted normal coor- 
dinates when r(no) = for points uq, for instance, belonging to a line on 
U. Such conditions can be satisfied for decompositions of metrics and con- 
nections on a small parameter like it is explained in details in Ref. [9j (see 
decompositions on a small eccentricity parameter e in Section [5]). Other 
possibilities can be found if we impose nonholonomic constraints, for in- 
stance, of type hi = but for nonzero /i4(x^,y^) and/or hl*{x^,y^); such 
constraints can be solved in non-explicit form and define a corresponding 
subclass of N-adapted frames. Considering further nonholonomic deforma- 
tions with a general decoupling with respect to a "convenient" system of 
reference/coordinates, we can deform the equations and solutions to config- 
urations when terms of type F • F became important. 

Theorem 2.3 ("Non quadratic" decoupling) The Einstein equations 
in GR (for instance, in the form (A.l^) and liA.13\) ). via nonholonomic 



9 



frame deformations to a metric g Iil0\) and matter source M. j5|) . when 

contributions from terms of type T ■ T are considered small for an open 
region [/ C V, can be transformed equivalently into a system of PDE with 
h-v-decoupling: 

Bl S2 ~1 r g'g* (.9*) , „// .91^2 (g'l)^ 1 

ni= H2 = [92 - -z ^ h gi - — J = - 1 , (24) 

25152 251 2g2 2g2 2gi 



' 2h3hi^^ 2/14 2h3^ 2/13/14^-3 2/i3 2h^^ 

= -T-T, (25) 

g _ Wfc (fel)^ ^5^4 1 , K ( dkhs ^ dkhA dkhX 

+ 2/1,^'= +U4-' 2-3j2/,4-°' 

p ^ :^fc ri^oo (^3)^ Ihhl hi f dkhj, dkhA dkh^ 

2/13 2/13 2/14 ^ 4/13^/13 /i4 y 2/13 

= (9i - Wi) In |/i4|, (9fe - Wk)wi ^ (di ~ Wi)wk,n* = 0, dink = 9^^^, (28) 



and 



w° = 0, drWk = dkWi,TT° (a^ - In |/i3|, (Sfe - nfe)n, = (9, - njrifc, 
efccj = Sfeo; — {wi + Wj) w* — {ni + tiJ w° = 0. (29) 



Proof. It is a consequence of Conclusion 12.11 and Theorems 1 2 . 1 1 and 
for superpositions of ansatz (f22l) and (j23]) resulting into ([9]). If we repeat 
the computations from Appendix [B] for geometric data (jlOp considering that 
contributions of type V ■ V are small, we can see that the equations (f25]l ~ (f27l) 
are derived to be respectively equivalent to sums of (fT2]) - (fH|) and (fT7|) ~ (fT9|) . 
The torsionless conditions (|28p consist a sum of similar conditions psp and 

dsni). □ 

In general, the solutions defined by a system (p^ - ([29|) can not be trans- 
formed into solutions parametrized by an ansatz (j22p and/or (]23p . As we 
shall prove in Section [3l the general solutions of the such systems of PDE are 
determined by corresponding sets of generating and integration functions. 

A solution for (I24p ~ (l29p contains a larger set of /i ?;-generating functions 

than those with some N-coefficients stated to be zero. 



10 



2.2 Splitting of Einstein— Yang— Mills— Higgs equations 

In terms of the canonical d-connection D such nonholonomic interactions 
are described by PDE, 

87rG(^T;35+ (30) 

^ie[^,D''<^], (31) 
A( cDfgj - $2)$, (32) 

where the source is determined by the stress-energy tensor 

""Tps = Tr[^ (i?5$ D0<i> + Di3'i> Ds<P) - ^gp5D^<S? - g^iV($),(33) 

^""'Tps = 2Tr (^g^'-Fp.Fs, - ^gpsF.^F^^^ . (34) 

These equations can be derived foUowing a variational principle for a grav- 
itating non-Abelian SU(2) gauge field A = A^e^ coupled to a triplet Higgs 
field $ as in [151 IS]- For our purposes, the operator D is used instead of 
V and all computations are performed with respect N-adapted bases ([3]) 
and The nonholonomic interactions of matter fields and constants in 
([30]) ~(f32]) are treated as follows: The gauge field with derivative Z)^ = 
+ie[A^,] in a covariant gravitational background should be changed into 
Vs = T>s + ie[As,]. The curvature of vector field A^ is 

= e/jA^ - e^A^ + ie[Afs,Afj], (35) 

where e is the coupling constant, = —1, and [•, •] is used for the commuta- 
tor. The value $[o] in (p2]l is the vacuum expectation of the Higgs field which 
determines the mass ^ M = ^/Xt], when A is the constant of scalar field 
self-interaction with potential V{^) = ;jATr(<I>^Qj — <^^)^; the gravitational 

constant G defines the Plank mass Mpi = 1/VG and it is also the mass of 
gauge boson, ^ M = ev. 

For a series of assumptions on nonholonomic constraints, we can decou- 
ple the system of PDE (|30 p - (|32p using corresponding splitting for Einstein 
equations considered in Theorems 12.11 12.21 or I2.3[ 

Condition 2.1 To construct new classes of solutions of off-diagonal EYMH 
equations we consider that a "prime" solution is known for the system 
^3^-^3^ (given by data for a diagonal d-metric °g =[ °gi{x^), °ha{x^), 
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= 0] and matter fields °A^{x^) and °^{x^), for instance, of type con- 
structed in Ref. 118^ (see also Appendix \B.5\) ). We suppose that there are 
satisfied the following conditions: 

1. The "target" d-metric with nontrivial N- coefficients, for °g — t- ''g 
is parametrized in a form 

g - v^{x'')°g,{x^)dx'(S>dx' + 'na{x'',y'')°ha{x\x^')e''^e'' 

= g,{x^)dx' ®dx' +uj^{x^,y^)ha{x'',y'')e'' ®e'', (36) 
= dy^ + [wi + wi\dx\ e'' = dy^ + [ui + n^dx\ 



2. The gauge fields are nonholonomically deformed as 

A^{x\y^)= °A^{x^)+ ^A^{x\y'^), (37) 

where °^^(x^) is defined by formula jlB. 8\} and '^^^(a;*, y") are any 
functions for which 

F^^ = °Fp^{x') + '^Fp^{x\y'') = s^\Ep^, (38) 

for s = const and epfj, being the absolute antisymmetric tensor. The 
gauge field curvatures Fp^j^, °Fp^ and ^Fp^j, are computed by intro- 
ducing W.9(l and (3l^ into [3^). An antisymmetric tensor Fg^ (33\) 
satisfies the condition D^{\/\g\F^^^) = 0, which always give us the 
possibility to determine '^Fp^, ^A^, for any given °A^, °Fi3^- 

3. The scalar field is nonholonomically modified °^{x^) — t- <I>(x*,y") = 
^ri{x^,y°') °(^[x^) by a polarization ^rj is such way that 

Df,<P = and Hx^y") = ±«>[o]. (39) 

This nonholonomic configuration of the nonlinear scalar field is not 
trivial even with respect to N-adapted frames V($) = and ^Tps = 0, 
see formula Ii33\) . For ansatz Ii36\) . the equations Ii39\) transform into 



{d/dx'-Ai)^ = {wi + Wi)^* + {ni + ni)<^°, (40) 
{d/dy^-A3)<^ = 0, (9/9y^ - A4) $ = 0. 

So, a nonolonomically constrained/ deformed Higgs $ field (depending 
in non-explicit form on two variables because of constraint i39\) ) mod- 
ifies indirectly the off-diagonal components of the metric via Wi + 
and Hi + and conditions (JU^ for ^A^. 
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4- The gauge fields i38\) with the potential ( |37[ j modified nonholo- 
nomically by <I> subjected to the conditions i39\} determine exact so- 
lutions of the system [25\) and (2B\) if the spacetime metric is cho- 
sen to be in the form i36\) . The energy-momentum tensor is com- 
pute^ yMjTQ _ _4g2^a_ Interacting gauge and Higgs fields, with 
respect to N~adapted frames, result in an effective cosmological con- 
stant = Svrs^ which should be added to a respective source ^A.15\) . 

Now we can formulate a Claim for decoupling PDEs for generic off- 
diagonal EYMH systems: 

Claim 2.1 An ansatz g = [qi °gi,ria °ha\Wi,ni\ \3b]) and certain gauge- 
scalar configurations (A, <I>) subjected to Conditions \2. 1\ define a decoupling 
of the EYMH system ^3C^) - [3^I in a form stated respectively by Theorems 
\2.1l \2.Sl or \2.3\ if the sources IIA.15\) are changed in the form 

= diag[r^] ^ + ^^'T^^ = diag[r^ - 4s'6'^]. 

Such a statement is a straightforward consequence of above assumptions 
1-4 when in N-adapted frames the contributions of matter fields is defined 
by an effective cosmological constant In a particular case, we can state 
Yq, = if there are not other energy-momentum sources. This statement 
on nonholonomic gravitational and matter field interactions is formulated as 
a Claim and not as a Theorem. This follows from our experience that not 
all EYMH systems can be modelled geometrically as nonholonomic Einstein 
spaces with effective cosmological constant. The solutions of such PDE are 
defined by very general classes of off-diagonal metrics but the matter filed 
equations and their solutions are substantially constrained nonholonomically 
and N-adapted to corresponding frames of references in order to encode the 
data into a nontrivial effective vacuum gravity structure. 

The main conclusion of this section is that generic off-diagonal Einstein 
spaces can be generalized to EYMH configurations with effective cosmolog- 
ical constant. Such generalizations preserve the decoupling property of the 
Einstein equations and impose certain classes of nonholonomic constraints 
and additional (Pfaff type) first order systems of PDE for matter fields. 
With respect to coordinate frames, such systems of equations describe very 
complex, nonlinearly coupled gravitational and gauge-scalar interactions. 

^such a calculus in coordinate frames is provided in sections 3.2 and 6.51 in Ref. [19| 
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3 Generic OfF-Diagonal Solutions for EYMH Eqs 



The goal of this section is to show how the decouphng property of the 
Einstein equations allows us to integrate such PDE in very general forms 
depending on properties of coefficients of ansatz for metrics. Some similar 
theorems for 4-d, 5-d and higher dimension modifications of GR where 
proven in Refs \10\ [TT\ [9]. In this work, we generalize those results for 
EYHM systems following Corollary 12.11 and Claim 12. 1[ 

3.1 Generating solutions with weak one Killing symmetry 

We prove that the Einstein equations encoding gravitational and YMH 
interactions and satisfying the conditions of Theorem 12.11 can be integrated 
in general forms for h* ^ and certain special cases with zero and non-zero 
sources (jA.lSp . In general, such generic off-diagonal metrics are determined 
by generating functions depending on three/four coordinates. The bulk of 
known exact solutions with diagonalizable metrics and coefficients depending 
on two coordinates (in certain special frames of references) can be included 
as special cases for more general nonholonomic configurations. 

3.1.1 (Non) vacuum metrics with h* ^ 

For ansatz ([9]) with data = l,h^ = l,Wj^ = and = for (jlOp . when 
h* / 0, and the condition that the source 

T^, = diagir^ : T\ = r\ = Tix" ,y^) - As^;r\ = r\ = ^T{x'')-4s% 

(41) 

is not zero, the solutions of Einstein eqs can be constructed following 

Theorem 3.1 The EYHM equations of type ill]) - |-?-^[ ) with source |^ j[ ) 
can be integrated in general forms by metrics 

g = eie^(^')(ix^®dx^ + /^3(x^y=^)e3®e3 + /l4(x^y^)/^4(2;^/)e^®e^ 
e3 = dy^ + Wi(x^,y^)dx\ e^ = dy^ + ni{x^)dx\ (42) 

with coefficients determined by generating functions il){x''),(j){x^ ,y^), cf)* ^ 
Q,ni{x^) and h^{x'^,y^), and integration functions ^4>{x^) following recur- 
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rent formulas and conditions 



\T-4s 



Wi 



"4(A-4s2) 

|/l4| 

4(A-4s2)^ 

dic, 



(43) 

dy^, or (44) 
, i/ T = A = const ^ 4s^; 



2|T-4s2 



(ln|/i4|)* or 



ifT = A = const / 4s^ 



where constraints 



{di - Wi) \n\hi\,dkWi = diWk, 
dkh^,dink = dkUi. 



(45) 
(46) 

(47) 



must he imposed in order to satisfy the zero torsion conditions M5]) : we 
should take respective values = ±1 and it in |^^[ ) and/ or |^5[ j if we wont 
to fix a necessary spacetime signature. 



if /14 = 1. 



Proof. We sketch a proof which transforms into similar ones in (TUIIIIIE] 

• A horizontal metric gi{x^) is for 2-d and can be always represented 
in a conformally flat form ejc'^^^'^^dx* ® dx^ . For such a h-metric, the 
equation ([TT]) is a 2-d Laplace/wave equation (|13]) which can be solved 
exactly if a source ^T(x*^) — 4s^ is prescribed from certain physical 
conditions. 

• If /i| 7^ 0, we can define nontrivial functions 



In 



7 := 



, ai = hldicP, p = hl4>* (48) 



for a function (f)[x^ ,y^). If 0* 7^ 0, we can write respectively the equa- 
tions (|12p - (|14|) in the forms, 



(t)*hl 



.4 - 2/13/14 [T - 45^ 
I3wi + Oj = 0, 



(49) 
(50) 



^i +ini 



0. 
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For the last equation, we must take any trivial solution given by func- 
tions ni{x^) satisfying the conditions diUj = djUi in order to solve 
the constraints p^ . Using coefficients (|18]) with / and /3 7^ 0, 
we can always express Wi via derivatives of cj), i.e. in the form (j46p . 
We can chose any y^) with (/>* 7^ as a generating function and 
express /i4 and, after that, /13 as some integrals/derivatives of func- 
tions depending on and source T{x^,y^) — 4s^, see corresponding 
formulas (j44p and ()45p . The integrals can be computed in a general 
explicit form if T{x^,y^) = A = const 7^ 4s^, when other possible 
matter field (additionally to the considered YMH ones) interactions 
are approximated by an energy-momentum tensor as a cosmological 
constant, 

T^, = S^,T = 5^,A. (51) 

• A possible dependence on is present in function which must 
satisfy conditions of type ()B.6P in order to be compatible with (jlSp . It 
is not possible to write in explicit form the solutions for the zero torsion 
condition if the source T''^ is parametrized by arbitrary functions. 
Nevertheless, if X'^^ is of type ([5T|) . we get /14 ~ e^"^ and Wi ~ di4)/(j)* 
positively solve the constraints 

w* = [di — Wi) In I/14I and diWj = djWi, (52) 

transformed into (j)** di(j) — (p* {di(j))* = 0. By straightforward computa- 
tion, we can check that (fT5]) are satisfied by ([T7|) when n* = and Wi 
is determined by (f46]) . □ 



The solutions constructed in Theorem 13.11 and those which can be de- 
rived following Corollary 12.11 are very general ones and contain as particular 
cases (pehaps) all known exact solutions for (non) holonomic Einstein spaces 
with Killing symmetries. They also can be generalized to include arbitrary 
finite sets of parameters as it is proven in Ref . |9] . 

Remark 3.1 ( Example) Introducing H = — ln|/i4|, the sistem of equa- 
tions (4^ ! iSW ^'^'^ ^SW ca?^ be written in the form 

hs = -r F72T2, (53) 
(j)*w^ + d^(f> = 0, (54) 
w* + diH-WiH* = 0. (55) 
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Considering H{x^,y^) as a generating function, we can integrate 1155]) in 
explicit form using parametrizations wi =i Wi{x^ ,y^)x 2Wi{x^ ,y^) (in this 
formula, we do not consider summation on i). The function cj){x^,y^) can 
he any one for which the system of first order PDF ( [5x^| j has nontrivial 
solutions for any found w-i. So, we can define for any prescribed source 
T2 in /f53|). We conclude that the LC-conditions i flJ)) (relating Wi to /14, 
which does not allow a "complete" decoupling) can he solved also in very 
general forms by further fixing of parametrizations, classes of functions and 
boundary conditions for H, 4> and T2. An explicit example of exact solutions 
can be derived by taking derivative on * for jj^p and using a subclass of 
functions (j) when {di(j))* = di{(j)*). We obtain Wj^j^ + w* + = 0. This 
system is equivalent/compatible to i f55)) if, for instance, for <I> := ln|(/)*|, 

H* = -^*, diH = di^. 

Such equations are satisfied by any functions H{ ~q) and $( "'"<;^), where 
:= ±v — x^ — x^ . 

3.1.2 Effective vacuum EYMH configurations 

We can consider a subclass of generic off-diagonal EYMH interactions 
which can be encoded as effective Einstein manifolds with nontrivial cos- 
mological constant A = 4s^. In general, such classes of solutions depend 
parametrically on A — and do not have a smooth limit from non- vacuum 
to vacuum models. 

Corollary 3.1 The effective vacuum solutions for the EYHM equations 
with ansatz for metrics of type j^^p with vanishing source are paramet- 
rized in the form 

g = e^e'''(^')(i2;^®dx^ + /^3(x^y3)e3®e3 + /^4(x^y3)/^^(3;'=,y4)e4 0e^ 
e^ = dy^ + Wi{x^,y^)dx\ e^ = dy^ + ni{x^)dx\ (56) 

where coefficients are defined by solutions of the system 

t/i + V" = 0, (57) 
(t>* hi = 0, (58) 
pwi + ai = 0, (59) 

where coefficients are computed following formulas for nonzero (p* and 
h\ with possible further zero limit; such coefficients and /13 and h^ are sub- 
jected additionally to the zero-torsion conditions ( [^7| j. 
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Proof. Considering the system of equations (jlip - (jl4p with zero right 
sides, we obtain respectively the equations ([ST]) - ([59|) . For positive signa- 
tures on h-subspace and equation ([57|) . we can take ^ = 0, or consider a 
trivial 2-d wave equation if one of coordinates is timelike. There are two 
possibilities to satisfy the condition The first one is to consider that 

/14 = h4{x^), i.e. /i| which states that the equation (j58p has solutions with 
zero source for arbitrary function hs{x^,y^) and arbitrary N~coefficients 
Wi{x^,y^) as follows from (j48p . For such vacuum configurations, the func- 
tions and can be taken as generation ones which should be constrained 
only by the conditions (jl7|) . Equations of type ([52]) constrain substantially 
the class of admissible Wi if hi depends only on x^. Nevertheless, can be 
an arbitrary one generating solutions which can be extended for nontrivial 
sources T and systems i^- 1^ and/or 1^- (|29D . 

A difi^erent class of solutions can be generated if we state, after corre- 

- 0, 



hl/^/\h3h4 



const. 



sponding coordinate transforms, (j) = In 

0. For such configurations, we can consider /i| 7^ 0, and solve (j58p as 

vW= °/i(v^)*, (60) 

for = const ^ 0. Such v-metrics are generated by any /(x*,y^), /* 7^ 0, 
when 

hA = -f {x\y^) and = ( %f [/* {x\y^)]\ (61) 

where the sines are such way fixed that for Nf — )• we obtain diagonal 
metrics with signature (-I-, -|-, -|-, — ). The coefficients = /3 = in (f59]l and 
Wi{x^,y^) can be any functions solving ()17|) . This is equivalent to 

w* = 29,111 I/I -2u;,(ln I/I)*, (62) 
dkWi-diWk = 2{wkdi-Widk)la\f\, 

for any ni{x^) when diUk = d^ni. Constraints of type nkh^ = dkh^ ()B.6P 
have to be imposed for a nontrivial multiple h^. □ 

Using Corollary 12. H the ansatz (f56]l can be "dualized" to generate ef- 
fective vacuum solutions with weak Killing symmetry on d/dy^. Finally, we 
note that the signature of the generic off-diagonal metrics generated in this 
subsection depend on the fact which coordinate x^,x^,y^ or is chosen to 
be a timelike one. 

3.2 Non-Killing EYMH configurations 

The Theorem 12. 21 can be applied for constructing non-vacuum and effec- 
tive vacuum solutions of the EYMH equations depending on all coordinates 
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without explicit Killing symmetries. 



3.2.1 Non vacuum off diagonal solutions 

We can generate such YMH Einstein manifolds following 
Corollary 3.2 An ansatz of type \22^) with d-metric 

g = ue^d.^^d.^+.^[±^^^^^e^^e^±^^^±^^ 

= dy^ - {d^(f>/<p*)dx\ e'^ = dy'^ + n,{x^)dx\ 

where the coefficients are subjected to conditions |?3[)~(f7^ o,f^d d^uj + 
{di<j)/(j)*)uj* — niUj° = 0, defines solutions of the Einstein equations R^p = 
(A — 4:s'^)ga/3 with nonholonomic interactions of YMH fields encoded effec- 
tively into the vacuum structure of GR with nontrivial cosmological constant, 
A - 4sV 0. 

Proof. It is a consequence of Theorems 12.21 and Corolarry 13.11 □ 

In a similar form, we can generate solutions of type (j23p when the confor- 
mal factor is a solution of d^uj — ui^{x^)uj* + {di4> / = with respective 
"dual" generating functions uo and (/> when the data (I43p - ()47p are re-defined 
for solutions with weak Killing symmetry on d/dy^. 

3.2.2 Effective vacuum off diagonal solutions 

Vacuum Einstein spaces encoding nonholonomic interactions of EYMH 
fields can be constructed using 

Corollary 3.3 An ansatz of type \2S^) with d-metric 

g = eie^f'^'^'^Ux' ®dx' + uj'^{x'',y'')[{^hf [f* [x'^y"^)]^ (^e^^ 

{x\y'')Ux\y^)e''®e% 
e^ = dy'^ + Wi{x^ ,y'^)dx\ = dy^ + ni{x'^)dx\ 

where the coefficients are subjected to conditions fgOj] -(fgl| ] . fTTj ) and dkU — 
Wild* — niUJ° = 0, define generic off-diagonal solutions of Rap = 0. 

Proof. It is a consequence of Theorem 12.21 and Corollary 13.11 □ 

Solutions of type (j23|) can be defined if the conformal factor is a solution 
of dkU) — w^{x^)u;* + n^{x^ ,y'^)uj° = with respective "dual" generating 
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functions uj{x^,y°') and (f){x^,y'^) when the data (|60p - (|62p and ()47p are re- 
defined for ansatz with weak Killing symmetry on d/dy^. 
Summarizing the results of this section, we formulate 

Claim 3.1 All generic off-diagonal non-vacuum and effective vacuum so- 
lutions of the EYMH equations determined respectively by Theorem \3.1\ and 
Corollaries \ 3.1{ \3.S\ and \3.3\ can be generalized to metrics of type For 
such constructions, there are used nonlinear superposition of metrics and 
their "duals" on v-coordinates in order to define non-Killing solutions of 
respective systems of PDEs from Theorems \2.S\ and/or [273\ 

The statements of this Claim are formulated following our experience on 
constructing generic off-diagonal vacuum and non-vacuum solutions. For 
such nonlinear systems, it is not possible to formulate ceratin general unique- 
ness and exhaustive criteria. Almost sure, not all solutions of EYMH equa- 
tions can be constructed in such forms, or related to any sets of prescribed 
solutions via " nondegernerate" nonholonomic deformations. The length of 
this paper does not allow us to present all technical details and general for- 
mulas for coefficients for ansatz for d-metric^ which can be constructed 
following this Claim. Explicit examples supporting our approach are given 
in sections [5] and [6l 

4 The Cauchy Problem and Decoupling of EYMH 
Equations 

The gravitational interactions of EYMH systems studied in this work 
are described by off-diagonal solutions of 

Ro,p = {A-As^)gap, (63) 

which can be found in very general forms with respect to N-adapted frames 
for certain nonintegrable spacetime 2 + 2 splitting of type N ([2]). An ef- 
fective cosmological constant (A — 4s^) encodes a gravitational "vacuum" 
cosmological constant A in GR and is induced by nonholonomic dynam- 
ics of YMH fields. Pseudo-Riemanian manifolds with metrics ga/^ adapted 
to chosen non-integrable distribution with 2 + 2 splitting and satisfying ()63p 
are called nonholonomic Einstein manifolds. Hereafter we shall refer to such 

*such formulas are, for instance, of type (|43ll- (l47p . with functions /i3(., j/^) and /i4(., y^) 
and further "duaHzation" /13 Z14 and /14 h-i{-,y'^) with corresponding re- 
definition of N-connection coefficients 
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systems of PDEs as nonholnomic vacuum spacetimes, regardless of wether 
or not an (effective) cosmological constant vanishes or can be "polarized" by 
gravitational and/or matter field interactions into some N-adapted diagonal 
sources admitting formal integration of gravitational field equations. 

The equations , and their N-adapted equivalents (|A.12p - (jA.13p , con- 
stitute a second-order system quasi-linear PDE for the coefficients of space- 
time metric g = {go/?}. This mean that given a manifold V of necessary 
smooth class such a quasi-linear system is linear in the second derivatives 
of the metric and quadratic in the first derivatives d^gaj3 (the coefficients of 
such PDE are rational functions of Qap)- To be able to decouple and formally 
integrate such systems is necessary to consider special classes of nonholo- 
nomic frames and constraints. For GR, this type of equations do not fall in 
any of the standard cases of hyperbolic, parabolic, or elliptic systems which 
typically lead to unique solutions. It is important to formulate criteria when 
such general solutions would be unique ones with a topology and differential 
structure determined by some initial data. How the differomorphysm, or co- 
ordinate, invariance and arbitrary frame transforms (principle of relativity) 
would be taken into account for N-splitting? 

In mathematical relativity [Il[3llH[5l[6l[71[8], it was proven a fundamental 
result (due to Choquet-Bruhat, 1952) that there exists a set of hyperbolic 
equations underlying ([63j) . The goal of this section is to study the evolution 
(Cauchy) problem for the system ()A.12p - (lA.13p in N-adapted form and 
preserving the decoupling property. 

4.1 The local N— adapted evolution problem 

In the evolutionary approach, the topology of spacetime manifold is cho- 
sen in the form V = M x ^V, where is a. 3-d manifold carrying initial 
data. It should be noted here that there is no a priory known natural 
time-coordinate even we may fix a signature for metric and chose certain 
coordinates to be time or space like ones. 

Definition 4.1 A set of coordinates {u^ = (x*,y")} is canonically N-har- 
monic, i.e. it both harmonic and adapted to a splitting N (0), if each of the 
functions satisfies the wave equation 

□n'^ = 0, (64) 

where the canonical d'Alambert operator □ := acts on a scalar 
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f{x, y) in the form 

for a d-metric g^/j = {gij,gab) (Q^ defined with respect to N-adapted (co) 
frames (Oj-j^j and canonical d-connection Dq, IIA.3\} . 

We can say that four such coordinates = {x^,y^) are N-adapted 
wave-coordinates. 

Lemma 4.1 In canonical N-harmonic coordinates, the Einstein equations 
re-defined in canonical d-connection variables iA.12\) can be written 
equivalently 



-2(A - V)g"'^ = 0, (65) 

i.e. such PDE for g'^^ (using algebraic transforms, for ga/s) form a system 
of second-order quasi-linear N-adapted wave-type equations. 

Proof. It is a standard computation with respect to N-adapted frames 
by using formulas ()A.3p . ()A.9P and ()A.7p . If the zero torsion conditions 
(IA.13P are imposed, we get well known results from GR but (in our case) 
adapted to 2 + 2 non-integrable splitting. □ 

This Lemma allows us to apply the standard theory of hyperbolic PDE 
(see, for instance, |20J). Let us denote by Hj'^^ the Sobolev spaces of func- 
tions which are in L?'{K) for any compact set K when their distributional 
derivatives are considered up to an integer order k also in L'^{K). We shall 
also use N-adapted wave coordinates with additional formal 3 + 1 splitting, 
for instance, in a form = { tT), where is used for the timelike 
coordinate and are for 3 spacelike coordinates. "Hats" can be eliminated 
if such a splitting is considered for arbitrary local coordinates. Standard 
results from the theory of PDE give rise to this 

Theorem 4.1 The field equations (A.l^) for nonholonomic Einstein man- 
ifolds have a unique solution g"^ defined by PDE h65\) stated on an open 
neighborhood W C M x of O C {0} x with any initial data 

g'^^QM)^Hlf and ^-^^^^Htt\k>V2. (66) 
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The set U can be chosen in a form that (U, g°'^) is globally hyperbolic with 
Cauchy surface oH 

There is no reason that a solution constructed using the anholonomic 
deformation method as we considered in section [3] will satisfy the wave con- 
ditions (|64p even we may give certain initial data for an equation ()65p . In 
order to establish an hyperbolic (and evolutionary) form of the such non- 
holonomic vacuum gravitational equations we should re~define the Choquet- 
Bruhat problem (see details and references in [H El HJ [U [8]) with respect 
to N-adapted frames. Using a 3+1 splitting of N-adapted coordinates, 
we wirte g"^ = {g'^^,g^^) and = {et,ej) and consider the d-vector field 
n'^{x,y) of unit timelike normals to the hypersurface { = O}0 There is 
no loss of generality if we assume that on such a hypersurface we have 

g" = -1 and g*^ = 0. (67) 

We can state such conditions via additional N-adapted frame/coordinate 
transform for any d-metric ([7]) with prescribed signature. It is also possible 
to re-define the generating functions for any class of off-diagonal solutions 
with decoupling of Einstein equations in order to satisfy ([67|) . 

Another necessary condition for the vanishing of a Dn^ is that this value 
is stated zero at = 0. This allows us to express in N-adapted form the 
time like derivatives through N-elongated space type ones, 

et (^g^g*^) = -e^ ( V&g*^) (68) 

with N-elongated operators. So, the initial data from Theorem 14.11 can not 
be fixed in arbitrary form if we wont to establish a hyperbolic (evolutionary) 
character for nonholonomic vacuum Einstein equations, i.e. to satisfy both 
systems ([65]) and ([M]) . Really, the last system of fist order PDE allows us to 
compute the time-like derivatives dg^^(0,u^)/d *m |{ t2:=o}if gjj |{ t^=o} and 
dgjj/d *n 1 1 ts:=o}^^ave been defined. We conclude that the essential data for 
formulating a N-adapted evolution problem should be formulated for a 3-d 
space d-metric 

^^^g:=gjj{x,y)e'^e\ (69) 

^Choosing corresponding classes of nonholonomic distributions N ([2]), we can relax the 
conditions of differentiability as in Refs. 1 2X1 122) (we omit such constructions in this work). 

^°We do not use labels for coordinates like 0,1,2,3 because the decoupling property 
of the Einstein equations and general solutions can be proven for arbitrary signature, for 
instance ( — h ++) , (+ H h) etc and for any set of local coordinates u" with a = 1, 2, 3, 4. 
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where are N-elongated differentials of type ([3|), and its N-elongated ([3]) 
time-like derivatives. 

Using the Theorem 14.11 and above presented considerations, we get the 
proof of 



Theorem 4.2 // the initial data \66\) satisfy the conditions (61) and 
and the so-called Einstein N-adapted constraint equations, 



E«/3 + Aga/3 j = 0, (70) 

are computed for zero distortion in IIA.I]} . then the d-metric stated by Theo- 
rem \4.1\ defines solutions of the nonholonomic vacuum equations \6^) and/or 

This theorem gives us the possibility to state the Cauchy data for de- 
coupled EYMH systems and their generic off-diagonal solutions in order to 
generate N-adapted evolutions. 

4.2 On initial data sets and global nonholonomic evolution 

We adopt this convention for spacetime nonholonomic manifolds V = 
M X ^V, were is a N-adapted 3-d manifold, when a Whitney sum 
T = h © V is stated by a space like nonholnomic distribution '^N 
of type ([2]) and there is an embedding e : '^V — )• V. 

Using the d-metric g (I69p , we can define the second fundamental form 

K of a spacelike hypersurface ^^inV, K(X,Y) := g(Dxn, Y), VX £T^V. 

-1 



The unity d-vector n = n^e^ = n^e* = n^e* = (^-\/|(?**| j e*, with nor- 
malization g(n, n) = g"^nQ,n^ = g**(n()^ = —1, is time-future and normal 
to ^V. The value K = {Kjj} is the extrinsic canonical curvature d-tensor 
of ^V. Imposing the zero torsion conditions (lA.lSp . K — )• {Kjj}, where 

Kjj = -^g*° (ejg^-+ejg^j - eagjj^ nt are components computed in stan- 
dard form using the Levi-Civita connection, but with respect to N-adapted 
frames. We can invert the last formula and write dtg^j = 2(g^*'nt)~^ Kjj+ 
{terms determined by ga^ and their space-derivatives}. Such formulas show 
that Kjj and dtgjj are geometric counterparts on hypersurfaces { *u = 0}. 

For the canonical d-connection = (^Dtj'Djj and d-metric t^'g in- 
duced on a spacelike hypersurface in a Lorentzian nonholonomic manifold 
V, we can derive a N-adapted variant of Gauss-Codazzi equations 
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In these formulas, t^'i?*^ is the canonical curvature d-tensor of ['^Ig, 
is computed following formulas ()A.6P as the spacetime canonical d-curvature 
tensor, n is the timelike normal to hypersurface when the local N- 
adapted coordinate system is such way chosen that d-vectors ej are tangent 

to '^V. Contracting indices, introducing divergence operator div determined 
by N-elongated partial derivatives ([3]), trace operator tr and absolute differ- 
ential d, we derive from above equations and ()A.12p this system of equations: 

divK — d{trK) = Sir J, momentum constraint; (71) 

[3]^-2A- ^^^g^^K\k-^ + {trKf = levrp, Hamiltonian constraint; 

C(-F, f'^^g) = 0, Einstein constraint eqs; 

where ^^j-R is computed as the scalar (jA.SP but for '^^g. In a general context, 
we consider that is embedded in a nonholonomic spacetime V with 
induced data ^ ^y, [^]g,K, , we have J := — (n, •) and p := (n, n) . The 

term C{J-, t^lg) denotes the set of additional constraints resulting from the 
non-gravitational part, including nonholonomic distributions N ([2]). If in 
such a set there are included the zero torsion conditions (|A.13p . we can 
omit "hats" on geometric/physical objects if they are written in "not" N- 
adapted frames of reference. The equations (jTlh form an undetermined 
system of PDE. For 3-d, there are locally four equations for twelve unknown 
values give by the components of d-tensors '^^g and K. Using the conformal 
method with the Levi-Civita connection, see details and references in Ref. 
[23], we can study the existence and uniqueness of solutions to such systems. 
Above considerations motivate 

Definition 4.2 A canonical vacuum initial N-adapted data set is defined 
by a triple ^ ^V, [^lg,K^ , where ( [^]g,K) are defined as a solution of [71\ ). 

If the conditions ()A.13P are imposed, the data ( ['^lg,K) are equivalent 
to similar ones ( t'^'g, K) with K computed for the Levi-Civita connection. 
Covering "^V by coordinate neighborhoods 0„ of points u G ^V, we can 
use Theorem 14. II to construct globally hyperbolic N-adapted developments 
(J^u,Su) of an initial data set (Uu C ^V, [^]g,K^ as in above Definition. The 
d-metrics generated in such forms will coincide after performing suitable N- 
adapted frame/coordinated transforms on charts covering such a spacetime 
wherever such solutions are defined. We can patch all data {Uu, gu) together 
to a globally hyperbolic Lorenzian nonholonomic spacetime containing a 
Cauchy surface ^V. This provides a proof of 
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Theorem 4.3 Any N-adapted initial data ^ '^V, t'^lg,Kj of differentiability 

class H^^^ X H^, s > 3/2, admits a globally hyperbolic, N-adapted and unique 
development (in the sense of Theorem \4- 1\ and above considered assumptions 
and proofs). 

For the EYMH systems defined by assumptions in Condition 12.11 we 
encode the data on such gravitational and gauge-scalar interactions into 
the nonholonomic structure of certain effective Einstein manifolds. 



5 Anholonomic YMH Deformations of Black Holes 

Gauge-Higgs nonholonomic interactions can be parametrized in such 
forms that they define off-diagonal deformations (for instance, of rotoid 
type) of Schwarzschild black holes. In this section, we study such EYMH 
configurations when A + *A = 0. Nonholonomic deformations can be de- 
rived from any "prime" data ( °g, °A^, °<I>) stating, for instance a diag- 
onal cosmological monopole and non-Abelian black hole configuration in 
[18j . We can chose such a constant s for *A when the effective source is 
zero (if '^A < 0, this is possible for A > 0). The resulting nonholonomic 
matter field configurations = + ^A^ (f37|) . F^jy = S\/[g[e^jy (p8]) 
and $ = subjected to the conditions ()39p are encoded as vacuum 

off-diagonal polarizations into solutions of equations ([M|) ~ (j25]) . 



5.1 (Non) holonomic non— Abelian effective vacuum spaces 

This class of effective vacuum solutions are generated not just as a simple 
limit A+ '^A — 7- 0, for instance, for a class of solutions (I42p with coefficients 
(jl3|) - P5]) . We have to construct off-diagonal solutions of the Einstein equa- 
tions for the canonical d-connection taking the vacuum equations R^^ = 
and ansatz g (|36p with coefficients satisfying the conditions 



eiV"(r,0) + e2V (r,0) = 0; (72) 
/13 = ±e~^ '^~r — for a given /i4(r, 0, (^), (j){r,6,Lp)= ^(j) = const; 

f^i = Wi{r,9,ip), for any such functions if A = 0; 

^ni{r,e) + ^ni{r,e) J {hlf \h4\-^^^dv, if n* ^ 0; 
Wr,0), if n*=0. 



rii 



Effective vacuum solutions of the Einstein equations for the Levi-Civita 
connection, i.e of Ra^ = 0, are generated if we impose additional constraints 
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on coefficients of d-metric, for e ^ ''^ = 1, as solutions of (j47p . 

' , K^Q- (73) 
0; wl - w[ = 0, w* = 0; 

(74) 

The constructed class of vacuum solutions with coefficients subjected to 
conditions (f72]) - ([71]) is of type (f56]) for ([57|l ~ (f59]) . Such metrics consist a 
particular case of vacuum ansatz defined by Corollary 13.31 with /i4 = 1 and 
a; = 1. 

Here we note that former analytic and numeric programs (for instance, 
standard ones with Maple/ Mathematica) for constructing solutions in grav- 
ity theories can not be directly applied for alternative verifications of our 
solutions. Those approaches do not encode the nonholonomic constraints 
which we use for constructing integral varieties. Nevertheless, it is possi- 
ble to check in general analytic form, see all details summarized in Refs. 
[9| I101 I11] (and formulas from Appendix [B]), that the vacuum Einstein equa- 
tions ()24p - ()15p with zero effective sources, ''T — 4s^ = and T — 4s^ = 0, 
can be solved by above presented off-diagonal ansatz for metrics. 



wiW2 I In 



W2 



±4 



5.2 Non— Abelian deformations of the Schwarzschild metric 



We can consider a "prime" metric which, in general, is not a solution of 
Einstein equations. 



-d^^d^-r'^iC) di^^di}-r'^{^)sm'^i} dip^dip + -!jj'^{S,) dt^ dt. (75) 



We shall deform nonholonomically this metric into a "target" off-diagonal 
one which will be a solution of the vacuum Einstein equations. The non- 
trivial metric coefficients in (j75p are stated in the form 



gi = -1, 92 = 
for local coordinates 



dr 



1 



■d, = (p,y'^ = t, where 



(76) 



and tu'^{r) = 1 



2/io 



+ 



If we put e = with fiQ considered as a point mass, the metric (|75|) 
determines the Schwarzschild solution. For simplicity, we analyze only 
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the case of "pure" gravitational vacuum solutions, not considering a more 
general construction when e = can be related to the electric charge for 
the Reissner-Nordstrom metric. In our approach, e is a small parameter 
(eccentricity) defining a small deformation of a circle into an ellipse. 

We generate exact solutions of the system ([57|) ~ (j59]) with effective A + 
■^A = via nonholonomic deformations — )• ^g, when gi = rjigi and ha = 

rjaha and Wi,ni define a target metric 

= m{Od£, <E) + V2{0rH0 d-d + (77) 
r]3{£„-d,ip)r'^{£,)sm'^ 19 5ip ® 6ip - rii{S„'d,ip)w'^{C) 5t®6t, 
Sip = dip + wi{£„d,ip)d£, + W2{^,'d,(p)dd, 5t^dt + ni{£„d)d^ + n2(i,'d)d'&. 

The gravitational field equations for zero source relate the coefficients of 
the vertical metric and polarization functions, 

hs = hl{b* f = r]s{^,^,ipy{Os^^^ ^, h4 = -b'' = -m(.L^,^W(.0, (78) 



for {rj^l = (/io) 1^4/^3! [(\/F/4l)*] - these formulas, we have to chose /iq = 
const ( /iQ = 2 in order to satisfy the first condition (f7¥|) ). where ha are stated 
by the Schwarzschild solution for the chosen system of coordinates and 774 
can be any function satisfying the condition 7]^ 7^ 0. We generate a class of 
solutions for any function 6(^, t?, 99) with b* 7^ 0. For different purposes, it is 
more convenient to work directly with 7/4, for rj^ 7^ 0, and/or /i4, for /i| 7^ 0. 
The gravitational polarizations rji and r/2, when rji = r]2r^ = e^^^'^\ are 
found from (I24p with zero source, written in the form ip** + ip" = 0. 

Introducing the defined values of the coefficients in the ansatz (iTTl) . we 
find a class of exact off-diagonal vacuum solutions of the Einstein equations 
defining stationary nonholonomic deformations of the Schwarzschild metric. 



2 

TU^ (5(y5 (g) + rj^vj'^ St (g) St, 



= _e^((f^0(f^+ d-d (g) d-d) - 
Sip = dLp + wid£,+'W2d-d, St = dt + ^nid£,+ ^n2d'd. (79) 

The N-connection coefficients Wi{^,^,ip) and "'^nj(^, -i?) must satisfy the 
conditions (I74p in order to get vacuum metrics in GR. 

It should be emphasized here that, in general, the bulk of solutions from 
the set of target metrics do not define black holes and do not describe obvious 
physical situations. They preserve the singular character of the coefficient 
w'^ vanishing on the horizon of a Schwarzschild black hole if we take only 
smooth integration functions for some small deformation parameters e. 
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5.3 Linear parametric polarizations induced by YMH fields 

We may select some locally anisotropic configm'ations with possible phys- 
ical interpretation of gravitational vacuum configurations with spherical 
and/or rotoid (ellipsoid) symmetry if it is considered a generating function 

= q{i,'d,^)+EQ{i,d,^). (80) 

For simplicity, we restrict our analysis only with linear decompositions on a 
small parameter e, with < e << 1. 

Using dMI), we compute {h*f = [(v^)*]^ [1 + e-^L-(^/y^)*] and 



WHY 

the vertical coefficients of d-metric (j79p . i.e /13 and /14 (and corresponding 
polarizations 7^3 and 1^4), see formulas ([78]) 1^ We model rotoid configurations 
if we chose 

2ii{^,^,ip) qojr) 

q=l and g = ——r- sm.{ujQip + 990), (81) 

r A^JL'^ 

for ;u(^, ??, ip) = /xo+e//i(^, 9j) (supposing that the mass is locally anisotrop- 
ically polarized) with certain constants ^, uoq and (fQ and arbitrary func- 
tions/polarizations ^i(^, "i?, ip) and qQ{r) to be determined from some bound- 
ary conditions, with e being the eccentricityl^ This condition defines a 
small deformation of the Schwarzschild spherical horizon into an ellipsoidal 
one (rotoid configuration with eccentricity e). 

The resulting off'-diagonal solution with rotoid type symmetry is 

™*g = -e^{di(^di+ d'd®d'd) + {q + eQ) 6t0 6t 

-4 \i^/\q\r] ' [1 + e-^ie/^/\q\r] ® (82) 

S(p = dip + wid^ + W2d'd, 5t = dt + ^nid^ + ^n2d'&. 

The functions g(^, (p) and q{S,, ^, (p) are given by formulas ([HT]) and the 
N-connection coefficients Wi{^,'d,<p) and rij = ^ni{^,-d) are subjected to 
conditions of type ([7^ . 

tfit(;2 ( In | — I ) = W2—w'i, w*^Q] (83) 



W2 

01 wl-w[ = 0, < = 0; X(C,i?)- in^(C,i?) = 
and V'(C) being any function for which tp" + tp" = 0. 



^'^Nonholonomic deformations of the Schwarzschild solution (not depending on e) can 
be generated if we consider e = and b'^ — q and (&*)^ = ■ 

^■^We may treat e as an eccentricity imposing the condition that the coefficient hi — 
= ?;4(^, 1?, ip)-cu^{£,) becomes zero for data (|8H) if r+ ~ 2/io / ^1 + sm{ujoip + ¥'0) j • 



29 



6 Ellipsoidal EYMH Configurations and Solitons 



We can consider nonholonomic deformations for the EYMH systems for 
arbitrary signs of the cosmological constant A and an effective nontrivial 
source A + *A 7^ containing contributions of nonholonomic YMH config- 
urations. Such classes of solutions can be constructed in general form for 
a system (p^ - ([27|) and (jl7|) with coefficients of metric of type ([I2]) . Such 
metrics consist a particular case of non-vacuum ansatz defined by Corollary 
13.21 with = 1 and uj = 1. 

6.1 Nonholonomic rotoid deformations 

Using the anholonomic frame method, we can generate a class of solu- 
tions with nontrivial cosmological constant possessing different limits (for 
large radial distances and small nonholonomic deformations) than the vac- 
uum configurations considered in previous section. 

Let us consider a diagonal metric of type 

lg = dC'^dC + r^iC) d9^d0 + r^ (^) sin^ 6 dip ^ dip + (6 dt(S) dt, (84) 

where nontrivial metric coefficients are parametrized in the form 51 = 1, 52 = 
f'^iC)-, ^3 = r^(^)sin^t?, hi = \vj'^{^), for local coordinates = = 
'&,y'^ = ip,y'^ = t, with C = J dr/\q{r)\^ , and \w'^{r) = — o"^(r)g(r), for 
q{r) = 1 — 2m{r)/r — Ar^/3. In variables (r, 9, if) , the metric (j84p is equiv- 
alent to (fRTj) . 

The ansatz for such classes of solutions is chosen in the form 

= (Mi'<^'> (d^ ® + d0(gi d0) + hsi^, 0, ip) 5ip®5ip + hii^, e, ip) St ® St, 
Sip = dip + wi (^, 9, ip) d^ + W2 (C, 'p) d9, 
St = dt + ni (C, e, ip) d^ + n2 {£., e, p) de, 

for hs = -hUb*)^ = V3iC,0,ip)r^iOsmH, = = m{i,9,p,) y^w^^). 
The coefficients of this metric determine exact solutions if 



4>"{i,e) + 4^'{i,e) = 2{k+ ^a); 



(85) 



^3 




4 (A + «A) 



1 



<i>{m). 




if n* + 0; 
if n* = 0; 
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for any nonzero ha and /i* and (integrating) functions ^ni{^,9), "^riii^^^O)^ 
generating function (j){^,9,ip), and ^(j){S,,9) to be determined from certain 
boundary conditions for a fixed system of coordinates. 

For nonholonomic ellipsoid de Sitter configurations, we parametrize 



rot 
A 



g = -e^(«'^) {d^ ^dC+ dO® de) + {q + eg) St St 



{^/\q\r [l+e^^{e/J\q\r]S^(^ S^, 



-hi 

Sip = dip + wid^ + W2d'd, St = dt + nid£^ + n2d'd, 
where q = \ — ^ ni^'^'V) ^ ^ _ SoIl^ _)_ ^p^^^^ chosen to generate 

an anisotropic rotoid configuration for the smaller "horizon" (when /14 = 
0), r+ ~ 2 ^fi/ ^1 + £=^jT- sin(a;o<^ + y^o)j , for a corresponding (/^(r). 

We have to impose the condition that the coefficients of the above d- 
metric induce a zero torsion in order to generate solutions of the Einstein 
equations for the Levi-Civita connection. Using formula ()83p . for (p* 7^ 0, 
we obtain that (/>(r, ip, 6) = In |/i4/\/l^3^4ll ™ust be any function defined in 
non-explicit form from equation 2e^'^(j) = A + *A. The set of constraints for 
the N-connection coefficients is solved if the integration functions in (j85p are 

chosen in a form when W1W2 ^In |^|^ =10* — w'l for w* 7^ 0; — w'l = 

for w* = 0; and take rij = ^ni{x^) for ^n'i{x^) — ^n2{x^) = 0. 

In a particular case, in the limit e — t- 0, we get a subclass of solutions of 
type ()86p with spherical symmetry but with generic off-diagonal coefficients 
induced by the N-connection coefficients. This class of spacetimes depend 
on cosmological constants polarized nonholonomically by YMH fields. We 
can extract from such configurations the Schwarzschild solution if we select 
a set of functions with the properties </> — t- const, Wi ^ 0,ni ^ and 
/14 — )• w'^. In general, the parametric dependence on cosmological constants 
and for effective YMH contributions, in nonholonomic configurations, is not 
smooth. 



6.2 Effective vacuum solitonic configurations 

It is possible to construct off-diagonal vacuum spacetimes generating by 
3-d solitons as examples of generic off-diagonal solutions with nontrivial 
vertical conformal factor lo. We consider that there are satisfied the con- 
ditons of Corollary 13.31 with /14 = 1 for effective vacuum solutions (such 
configurations may encode EYMH fields) and the Cauchy problem is stated 
as in Section [3l 
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6.2.1 Solutions with solitonic factor 6t;(a;^,y^,t) 

We take u = r]{x^ ,t), when = t is a time like coordinate, as a 
solution of KdP equation 



± T]** + [dtr] + r]7]' + €7]'")' = 0, (87) 

with dispersion e and possible dependencies on a set of parameters 9. It is 
supposed that in the dispersionless limit e — )• the solutions are independent 
on and determined by Burgers' equation dtr] + r] r]* = 0. For such 3-d 
solitonic configurations, the conditions ()47p are written in the form 

eir/ = rj' + wi{x^ ,y'^)ri* + ni{x^)dtr] = 0. 

For rj' = 0, we can impose the condition W2 = and n2 = 0. 

Such vacuum solitonic metrics can be parametrized in the form 

g = e^^'^'^dx^ C^dx^ + dx'^ C^dx'^)+[r]{x\y^,t)]^ ha{x\y^) e^C^e"", 
= dy^ + wi{x'',y^)dx^, e'^ = dy^ + ni{x'')dx'^. 

This class of metrics does not have (in general) Killing symmetries but may 
possess symmetries determined by solitonic solutions of (|87j) . Alternatively, 
we can consider that r/ is a solution of any three dimensional solitonic and/ 
or other nonlinear wave equations; in a similar manner, we can generate 
solutions for uj = r/(x^,y^,t). 

6.2.2 Solitonic metrics with factor uj{x^,t) 

There are effective vacuum metrics when the solitonic dynamics does not 
depend on anisotropic coordinate y^. In this case uj = r]{x^,t) is a solution 
of KdP equation 

+ {dt^+v^ + e^"y = 0. (88) 

In the dispersionless limit e — )• the solutions are independent on x^ and 
determined by Burgers' equation dtf) + rj rj' = 0. 

This class of vacuum solitonic EYMH configurations is given by 

2g = e'^(^')(dxi®dx^ + dx2®dx2)+ ^^(x^^)j\a(x^y3) e'^^e", 
= dy^ + wi{x'',y^)dx'^, e'^ = dy'^ + ni{x'')dx^; 

the conditions (H71) are eirj = r/* + ni{x^)dtr] = 0, 02?? = rf + n2{x^)dtT] = 0. 
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It is possible to derive an infinite number of vacuum gravitational 2-d and 
3-d configurations characterized by corresponding solitonic hierarchies and 
bi-Hamilton structures, for instance, related to different KdP equations (j88p 
with possible mixtures with solutions for 2-d and 3-d sine-Gordon equations 
etc, see details in Ref. [25j. 

A Nonholonomic 2+2 Splitting of Lorentz Mani- 
folds 

In a general metric-affine manifold V is endowed with a metric 

structure g and an affine (linear) connection structure D (as a covariant 
differentiation operator). A linear connection gives us with the possibility 
to compute the directional derivative DxY of a vector field Y in the direction 
of X. It is characterized by three fundamental geometric objects: 

1. the torsion field is (by definition) T{X,Y) := Dy^Y - DyX - [X,Y]; 

2. the curvature field is n{X,Y) := DxDy - DyDyi - £'[x,Y]; 

3. the nonmetricity field is Q{X) := Dxg- 

Introducing X = and Y = e^, defined by ([3]), into above formulas, 
we compute the N-adapted coefficients and symmetries of = {T'^^^j and 
corresponding fundamental geometric objects, 

r rp^ / rpi rpi rpa rpa rpa \ \ , 

~ \^ al3 ~ \^ jk^^ ja^^ ji^ ^ bi^ ^ be) S i 
TZ = {R"i3.yS= {R\jk^R"'bjk^R\ja^R^bja^R\ba^R^bea}}'^ Q = {Q'^a/s}- 

Every (pseudo) Riemannian manifold {V, g) is naturally equipped with a 
Levi-Civita connection D = V = { ^1"^^^} completely defined by g = {gafs} 
if and only if there are satisfied the metric compatibility, ^ Q{X) = Vxg = 
0, and zero torsion, = 0, conditions. Hereafter, we shall write, for 

simplicity, ^T'^^^ = T^^i^- It should be emphasized that V does not preserve 
under parallelism and general frame/coordinate transforms a N-splitting ([2]). 
Nevertheless, it is possible to construct a unique distortion relation 

V = D + Z, (A.l) 

where both linear connections V and D (the second one can be considered 
as an auxiliary linear connection, which in literature is called the canonical 
distinguished connection; in brief, d-connection) and the distortion tensor 
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Z, i.e. all values in the above formula, are completely defined by g = {gap} 
for a prescribed N = {N^}, see details in [9l [TOl [TT]. 

Theorem A.l With respect to N-adapted frames ^ and the coeffi- 
cient of distortion relation II A . 1\) are computed 

r% = f% + z^„,, (A.2) 

where the canonical d-connection D = { T'^^^ = (L*^, L^^, Cj^, C^^)} is 
defined by coefficients 

L)k = ^5*'' (efc^ir + ejQkr - ergjk) , (A.3) 
Ltk = eb{N^) + ^g""' {ekgbc - gdc ebNi - gab ecNi 
Cjc = ^9'''ecgjk, Cj^c = 2^"'' (^<^9bd + ebgcd - edgbc) , 
and the distortion tensor „ is 



^jk 




Cjbgikg"'^ 


^bk 


_ 4 


-■^ab fj\c 
^cd ^kb': 






— "0.d, rpc 
^cb jd 


ih 

jk ~ 




- gikg'"") 



^'■jk-' ^bk — l^^'-jkdchg - ^jk ^hb-: 



2 i^' bk — 2 
1 

kb — 2^''jkScbg + "J 
^bc — ^-- ^ab 



(A.4) 



Tja.gcb + Tjf^gca 



2 

"^Ef^ = + gcdg"^). The nontrivial 

coefficients ^2"^ and T''^^^ are given, respectively, by formulas ^ and, see 
below, iA.5\) . 

Proof. It follows from a straightforward verification in N-adapted 
frames that the sums of coefficients (jA.Sp and ()A.4p result in the coeffi- 
cients of the Levi-Civita connection T'''^^ for a general metric parametrized 
as a d~metric g = [gij,gab] ©• □ 

All geometric constructions and physical theories derived for geomet- 
ric data (g,V) can be equivalently modeled by geometric data ^g, N,D 
because of unique distortion relation ()A.ip . 

Theorem A.2 The nonholonomically induced torsion T = {T''^^^} o/D is 
determined in a unique form by the metric compatibility condition, Dg = 0, 
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and zero horizontal and vertical torsion coefficients, T^j^ = and T^^^ = 0, 
but with nontrivial h~v- coefficients 

rpi ji ri rpi /^i rpa r\a rpc J- c. { AT^\ 

^ jk — ^jk ~ ^kjT^ ja — ^jb^ ^ ji — ji'-'-aj — ^ aj ~ ^aV^ j ) , be — '-'be ~ ^cb- 

(A.5) 

Proof. The coefficients (|A.5P are computed by introducing = D, 
with coefficients (|A.3p . and X = ea,Y = (for N-adapted frames Q) 
into standard formula for torsion, 7~(X, y) := DxY — DyX — [X,Y]. 

□ 

In a similar form, introducing D and X = eo,,Y = ep,Z = into 
TZ{X,Y) := DxDy — DyDx — -D[x,y]i we prove 



Theorem A. 3 The curvature TZ = {R-^g^^} of D is characterized by N- 
adapted coefficients 



R\jk 




- e]L\k - 


r ^ hj^ mk 


T m T i f^i r\a 
^ hk^ rnj ^ /la" fcjJ 




^ bjk 




- ej^'bfe ^ 


fc fa 
' ^ bj-^ ck 


re T a /^a qc 
^ bk^ cj ^ ^ 6c" fcj! 


(A.6) 


Si 

jka 


= eaL'jk 


- DkC\, 


1 y^i rpb 

+ ^ jb-^ kaJ 


R'^bka — ^aL'^bk ~ DkC\^ 


I + ^^d'^'^ka^ 


fii 

^ jbc 


= e,&jb 


- ebC\^ + 


■ jb^ he ~ 


'-^ jc'^ hbJ 




'na 

-fl bed 


= s.dC\^ 


~ ScC^bd 


r^e r^a 

bc^ ed ~ 


y^e /^a 
bd^ ec- 





We can re-define the differential geometry of a (pseudo) Riemannian 
space V in nonholonomic form in terms of geometric data (g,D) which is 
equivalent to the "standard" formulation with (g, V). 

Corollary A.l The Ricci tensor 'R.ap '■= '^"'a/^-y (HIliP o/D is character- 
ized by N-adapted coefficients 

— {Rij '■— R^ijk^ ^ia ■= ^R^ika^ ^ai '■— R^aib: ^ab R^abe}- {-^■'^) 

Proof. The formulas for /i-f-components (jA.7p are obtained by con- 
tracting respectively the coefficients ()A.6p . Using D (]A.3p . we express such 
formulas in terms of partial derivatives of coefficients of metric g ([T]) and 
any equivalent parametrization in the form d?]), or ([8]). □ 

The scalar curvature of D is by definition 

:= g'^^K^p = g'^Rij + g^'K^. (A.8) 
Using (|X7ll and (lASp . we can compute the Einstein tensor Eq,^ of D, 

= Ra/3 - ^ga/3 'R- (A.9) 
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In general, this tensor is different from that constructed using (|A.10p for the 
Levi-Civita connection V. 

Proposition A.l The N-adapted coefficients T'^^i^ o/D are identic to the 
coefficients T^'^^^ ofV , both sets being computed with respect to N-adapted 
frames ^ and if and only if there are satisfied the conditions L'^j = 
ea{Nf),C]^ = d and = 0. 

Proof. If the conditions of the Proposition, i.e. constraints ()A.13p . are 
satisfied, all N-adapted coefficients of the torsion T"*^^ (jA.Sp are zero. In 

such a case, the distortion tensor 'ZP^^ ()A.4p is also zero. Following formula 

()A.2p . we get T'*^^^ = P"^^^. Inversely, if the last equalities of coefficients are 
satisfied for a chosen splitting ([2]), we get trivial torsions and distortions of V. 
We emphasize that, in general, D 7^ V because such connections have differ- 
ent transformation rules under frame/coordinate transforms. Nevertheless, 
if r'*^^^ = P'^^^ in a N-adapted frame of reference, we get corresponding 
equalities for the Riemann and Ricci tensors etc. This means that the N- 
coefficients are such way fixed via frame transforms that the nonholonomic 
distribution became integrable even, in general, the frames ([3]) and (jH are 
nonholonomic (because not all anholonomy coefficients are not obligatory 
zero, for instance, = daN^ may be nontrivial, see formulas ([6])). □ 

In order to elaborate models of gravity theories for V and/or D, we have 
to consider the corresponding Ricci tensors, 

Ric={Rp^ := R%J, iorV = {T\i,}, (A.IO) 
and Ric = {R := for D = {f (A.ll) 

Theorem A. 4 The field equations in GR can be re-written equivalently 
using the canonical d-connection D, 

R. - = r^s, (A.12) 

Llj = CaiN'r), C% = 0, n^^^ = 0, (A. 13) 

where the scalar curvature '^R := g^*^!! and source d-tensor Yps is such 
way constructed that T^^ — )• xT^^ for D — )• V, where Tjsg is the energy- 
momentum tensor in GR with coupling gravitational constant x. 
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Proof. The Einstein equations are written in "standard" form for V, 
Eps = Rps - R = xTps, (A.14) 

where R := g^^R ps- On spacetimes with conventional h- and v-splitting, 
we can define geometrically (or following a N-adapted variational calculus 
with operators ([3|) and Q) a system of field equations with the Einstein 
tensor E^^ ()A.9p . written in the form ()A.12p . In general, both systems of 
PDE are different. But if the constrainst ()A.13P are imposed additionally on 
D, we satisfy the conditions of Proposition IA.lt when T'^^p = T"^^^ results 

in Rjis = R ^5 and = E^/^. The coefficients ()A.5P are computed by 
introducing D = D, with coefficients (lA.Sp . and X = eQ,,y = (see ([3])) 
into standard formula T{X, Y) := D^X - DyX - [X, Y]. □ 

We consider matter field sources in (IA.12P which can be diagonalized 
with respect to N-adapted frames, 

T^, = diag[r^ : Y\ - r\ - T(x^ y^) + T(a;^ y^); = r\ = -T{x'')]. 

(A.15) 

Such a formal diagonalization can be performed via corresponding frame/ 
coordinate transforms for very general distributions of matter fields. 

An very important result which can be obtained using the anholonomic 
deformation method [9l 110^ [TT]. and developed in this work, is that the Ein- 
stein equations for D (|A.12p decouple for parametrizations of d-metrics ([7]). 
The corresponding systems of PDE can be integrated in general forms with 
one Killing symmetry and for "non-Killing" configurations. This allows us 
to generate exact solutions of standard Einstein equations ()A.14p in GR 
imposing additional nonholonomic constraints ()A.13p . It should be empha- 
sized here that if we work from the very beginning with the Levi-Civita 
connection for metrics dHj) computed with respect to coordinate, or other 
not N-adapted, frames, we can not "see" a possibility of general decoupling 
and "formal" integration of the gravitational field equations. 

B Proof of Theorem EOl 

For oj = 1 and = const, such proofs can be obtained by straightfor- 
ward computations [9|. The approach was extended for u ^ 1 and higher 
dimensions in [10\ lllj. In this section, we sketch a proof for ansatz ([9]) with 
nontrivial depending on variable when u; = 1 in data (jlOp . At the 
next step, the formulas will be completed for nontrivial values w / 1. 



37 



If R\ = R\ and R\ = R\, the Einstein equations ()A.12p for D and 
data ()B.ip (see below) can be written for any source ()A.15p in the form 

E\ - E\ = -R\ = T(a;^y3) + Y(x^y^/), ^ e\ = -R\ = "T(a;'=). 

The geometric data for the conditions of Theorem 12.11 are g-i = gi{x^) and 

g^ = h^{x\y^),g^ = h^{x\y^)h^{x\y^),Nf = w,{x\y^),Nt ^n,{x\y^), 

(B.l) 

for /13 = 1 and local coordinates = = (x^, x^, y^, y^). For such 

values, we shall compute respectively the coefficients of in ([6]), canonical 

d-connection T^^^ ()A.3P , d-torsion T"'^^ (jA.SP , necessary coefficients of d- 
curvature R^q^^ (IA.6P with respective contractions for R^/s := '^aP'y ^^■'^1^ 
and resulting ^R (lA.Sp and E^,/q (IA.9p . Finally, we shall state the conditions 
(jA.lSP when general coefficients (jB.ip are considered for d-metrics. 

B.l CoefRcients of the canonical d— connection 

There are horizontal nontrivial coefficients of F'^^^ (jA.Sp , 

L)k = ^9'^ (efc.9ji + ejgki - eig^k) + ^.g'^ (efe5j2 + ej5fc2 - 625^7=) 
= ^5'^ (<9fc5ji + djgki - SigjTc) + ^5'^ (5fcgj2 + 9jgfc2 - d2gjk) , 
i-e- ijfc = ^ (c'fcSji + - g'k) ' 2jfc = ^ (^fe5'j2 + a^Sfcz - g'jk) ■ 

The h-V"Components L^f^ are computed following formulas 

Uk - e6(A^|) + -^-M;,3-<?3e6iV|-536(A^,')*-.g4fc(iV|r] = e6(iV3) 
2.93 

+ ^[afe5M - Nigt, - 7Vfc4g,°3 - 53 e.TV^ - .g3fc(iVfe')* - .946(^^^)1 

= eb{wk + Wk) + -^[dkgbs - {wk + Wfe).9b3 - (""-k + 2i/c)563 
-53 e6(«^fc +W/c) - gsbWk* ~ gAbnk*], 

Uk = e^b{Nt) + ^[ekgM - 34 ebNt - g^b^Nlf - ff4fc(iV,')°] = eb{Nt) 
254 

+ ^[dkgM - A^i^56*4 - A^fcffM - 54 e,7V4 - 53fc(iVi?)° - g^biNtf] 
^54 

= e6(nfc + Uf^) + — [cffeffM - {wk + Wfc)fffc4 ^ ("fc + Zife)3b4 
254 

-54 e6(nfe + nfe) - 53&uifc - 9ibn°k]- 
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In explicit form, we obtain these nontrivial values 



1 r *i h^h^ * * I -'^ r *i ^ 

2^3 2/13 2^4 2 



L,k - ^[9.54 - Wkg, - n.g^] = - - n. — . 

For the set of h~v C~coefficients, we get Cj^ = \g'^^'^§^ = 0. The v- 
components of C~coefficients are computed following formulas 

Cbc = ir~ (ecS-bs + ec5c3 - gL) > C^c = J- (ecff64 + e65c4 - 56c) > 
; „ r'S _ S3 _ '^s ris _ 93 _ n r'S _ sS _ ^1^4 ri4 _ S3 _ 

1- e. <--33 - - 2^:^' '-"34 - 2iJ - '-"44 " "2^^ " ^^''-"33 " "2^1 " 

0, C34 = ^ = ^,C'44 = ^ = Putting together the above formulas, 
we find all nontrivial coefficients, 

^11 — 7; — ! ^12 — 7^ — 1^22 — — J ^11 — t; — ' ^12 — 7; — J ''^22 — 7; — ; l^-^J 
2gi 251 2gi 2g2 2g2 2g2 

-4 _ dkjhih^) Wkhl _ X Zi4 f 3 _ ^fcfea w^/tg -3 _ ^4^4 * 

" 2/14/14 ^ ^"'^+-'=^2^'^3fe-^ ^'^^fe-Z^'^fc, 

f4 _ _ ^3 ;=i3 _ ^4^4 pi4 _ ^3^3 ^=,4 _ hi -4 _ Z^t 

^3^. - ^n,, C33 - — , C44 - — ^, C33 - C34 - — , C44 - — . 

We shall need also the values 



^3 — C^33 + C^4: — + Ci — Cf^ + C44 — (B.3) 



B.2 CoefRcients for torsion of D 

Using data (IRT]) for Wi = Ui = 0, the coefficients = ej (iVf ) -ej(iV°) 
([6]), are computed 

= dj (Tvf ) - ai(7v;) - Nf{N;^y - N^iN^^y + Nf{N^y + Nf{N^y 

= d^{Nf)-d,{N^)-w,{N^y + w,{Nty. 
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We obtain such nontrivial values 

17^2 = ^^21 = d2Wi — diW2 — W1W2 + W2W\ = w'l — w* — ^1^2* + W2w\; 
f2^2 = ^^21 = d2n\ — din2 — winj + W2n{ ~ n'^ — n* — WiTi^ + W2n{. (B.4) 

The nontrivial coefficients of d-torsion (|A.5p are T"^ = —^""ji (|B.4p and 
T^j- = L^j- — ea{Nj). We find zero values for other types of coefficients, 

rj-ii 'j^i 'ji r\ i-pi p^i r\ rpa x~Ya /^a ri 

jk ~ ^jk ~ ^kj — ^ ja — ^jb — ^ be — '-'be ~ ^ eb — ^■ 

We have such nontrivial N-adapted coefficients of d-torsion: 





- f3 




dkhs 

Wk 

2/13 


K 
2/13 


-Wk, 






- e4(iV|) - 


- 2h, 


^3k 


— f4 


tL 


— ^ik 




dk{hih^) 
2/14/14 


- Wk 


K 

nk- 

2h4 : 


Tf2 




- W* — W1W2 


+ W2Wl, - 


-T^2 


= n\ - n. 



hi 



(B.5) 



If all coefficients (|B.5P are zero, T'^^^ = T'*^^^ . 

B.3 Calculation of the Ricci tensor 

Let us compute the values Rij = R^jj. from (|A.lip using (|A.6p . 

^\jk = ^kL\j ~ ^j^hk + ^hjKnk ~ ^Tk^mj ~ ^hJ^'jk 

= dkL\j ~ 9jL\,^ + Uf^jU^j^ — L'^^.Ll^j, 
where we put = and 

ekLlj = dkL],^ + N^daLl^ = dk% + Wk {lIj)* + rik (K^Y = dk%. 
Taking derivatives of ()B.2p . we obtain 

f 1 _ ( 9* \' _ 9i' (g*) Q f 1 

o'l^ 11 - \7^) - n TT, — T2' i: 

251 2gi 2(91)'' 



51^2 = i-^r = -^ 



92 ^. _ 92 , .9l32 



251 251 2 (51)^ 



2<72' 2.92 2 (.92)' 



2.91 


5r 
2.91 


9'i9'i 

2(ffif 






ff'l 




9*i92 


252^ 


252 


2(32)' 


2^2 


2.92 


9292 
H92f 


7 
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For these values, there are only 2 nontrivial components, 



ri _ 92 dig' ia'f , 9i 5'i52 ia'if 

-n- 515 — "~ — ~o : r 



92 , , {9-2? 9'{ , 9W2 , (a'lf 



2^2 4^152 4(52)2 2^2 4(52)2 45152' 
Considering i?ii = —R'^112 and i?22 = -^^212; when 5* = l/5i, we find 

iii = i?i = -^[52"-^-^+5r-^-^], 

25152 251 252 252 251 



which can be found in equations (j25p . 

The next step is to derive the equations ([26ll . We consider the third 
formula in ()A.6p . 



bka 



QyO. V Q^k ^dk'-'ba ~ ^bk^da ~ ^ak^bd) "r ^bd-'-ka 

bk /~ic I /—(C rpd 
-^-^ - ba\k + ^ bd^ ka- 

bka dy'^ ba\k ^ ^bd^ka' 



Contracting indices, we get R,^ = R\^^ = ^ - Cl^^ + C^A- For 
Cb '■= Cha' '^^ write 

Cb\k = ^kCb - L%k^d = dkCf, — N^dgCb — L%k^d 
= dkCb — WkCl — UkCh ~ L\kCd- 

We split conventionally R^k = [i]Rbk + [2]Rbk + [3]Rbk, where 

[l]Rbk = (^Lff^j + (^Lfi^^ , [2]Rbk = —dkCb + WkCb + rikCb + L'^bk^d, 

ry /^a rpd Pf3 rp3 i Pi3 rpA , Pi4 rynS i r7n4 

[3\^bk — ^bd-l-ka — ^b'i-'-kd,^ '^bA^k'i + ^bd,'''kA + '^bA-'-kA- 
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Using formulas ()B.2p . (jB.SP and (|B.3p . we compute 



[1] ^3fc 



[2]R3k 



[3] ^3fc 



Ltkl = 



dkhs _ h3_ 
2h3 ""^ 2h3 



"■3 
2h3 



dkhs 



h* 



h* 

2h4' 



2h3 

dkh3 
2h, 



2/14 

3(M)! 

4(/i3)2 

dkhs 
2h, 



2h3 

hi* 
2hl 



h* h* 



"-4 > 
2fu' 



1 



, * ii4 



2{h4y 4/13/14^ 



1 /i^ ftr 



2 ''2/14 



^3 ^ii3 

/^3 rTi3 I y^3 rTi4 
^33-'A;3 + ^34-'fc3 



2/13 2/14 2 /13 /l4 



Pi4 ^3 
^33-'fe4 



' '^34-'- k4 



2/13 



^fc^3 

2fe 



2/13 

^3 c'fc/13 



- Wk 



' Wk 



2h3 

*\2 



+ 



Wk 



2/13 2/13 



2/14 V 2h 



4(/l4)2 
9fe/l4 , 5feft,4 



2hA 



rik 



dkijuk^ 
2/14/14 

/t4 Ikl 

2fH2h^ 



hi ft-4 



2h. 



Summarizing, we get 



Wk 



i{h: 



I? 



1 hi hi 



2/14 4 (/l4)2 



hi dkhs 



Idkhl 



4 /l3 /l4 
_^ 1 /l|9fc/l4 



^ Hh^^ ''2hi 2h^ 



2/14 2/13 2 /i4 4 (/i. 



^4 C^fcZk 

2/m 2/14 ■ 



which is equivalent to (j26p if the conditions n^/i^ = 5^/14, see below formula 
()B.6p . are satisfied. 

In a similar way, we compute R^k = [i]Rik+ [2]^4A; + [3]-R4fc; where 



[l]RAk 



f3 



+ [L\k 



[2] Rik 



-dkCi + WkCl + nkCl + P^kCs 



+L 4fcC4, [3]ii4A: — C^dTf 



a rpd 
ka 



^43-'- k3 "T '-^44-'A:3 "-"43-^/04 "r ^44''- k4- 



For the first term, we use respective coefficients of d-connection, L?, and 
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L^fc, from 



[l]R4:k 



f3 



** hi , 



^ik 



2(/l3)* 



1 oo 

(hi 



h* h° 
"-4 \ —4 \c 

dklki iLidkhji 



2h, 



In order to compute the second term, we use formulas (|B.3p . for C3 and 
C4, and p.2p . for L^^ and L|^, 



dk{hihi 



nk 



hi 



4«4 



"3 



2/14/14 



2/14 



"-fe 



/l4 



/l4 



'■ 2h3h^ 



2h^ J 2h 



3 

2/^3 



2/14 







/l4 


9fc/l4 /l4 


A2/14, 


/ 2/I4 


2^ 


2/14 2/I4 



"-4 
2/14 



9fc/i^ h° dkhl h°dkh 



2/14 2/14 



2/14 



2(/l4 



For the third term we use formulas (IB.2p . with C43, C44, CI3, CI4, and the 
formulas with f ^g, f^, f 3^, f,\, 



/-i3 'TIS I /-<3 '714 1 PiA /T-iS I y?~(4 ^J^A 
"-"43-^^3 "T '-"44-'fc3 "T '-"43-^fc4 + *-^44^fc4 



2/13 



1 



'2/14 V ^^"2/14 



2h 



+ 



9t/i, 



4 ^"-4 
(9fc/l4 /I4 



'fc«4 ^fc/ji4 
2/I4 



2/1, 



/14/1 



i±4 



2/14 



/I4 9fc/l4 

2^ 2/14 ■ 



2/14 2/14 

Taking the last equalities in above three formulas, we get 



2/13 



■nk) 



[l]R4k 



[2] Rik 



[3]^4fc 



-/I4 



Wk 



2/13 

"2/14 2/14 



"-4 

2/t 



/i|/i;5 



M /^:4 



3 2(/l3)2 
/l4/l4 / /I3 

2^^ 



= 2/13/1 



3^3 



, dkh4 hi dkhi hi 



2/14 2/14' 



2/14 2/1, 



2/1, 2/1, 



h4 



h* 
»4 

2/14 

9fc/i:4 
'2h7 



~2h^ ~ 2(/i4)2 ' 



2/l4 2/14/ ''4/l4^ 



^4 

2^ 



9/c/l4 /I4 

2/14 2^ ~ 2^ 



hldkhi 

° 9fc/l4 



it4 



2/I4 
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We summarize above three terms, 

R - -^n*f ^4 , KK Kh* hi\ 



-A- I —4 \ _ —4 —4 I / —4 

V 2/I4 2(/),4)2 V2/I4/ 2/I4 2/I4 ^2hJ J ' 2/I4 2(/i,4)2 

dklH^_ dkh^ hi ^ dkhl , kldkh^ _ 5fcfe4 ^4 _ ^_dkh4_ 
2hi 2h^ 2h^2h^~ 2h^ 2{h^Y ~ 2h4 2h^ ^ 2h^ 2h^ ' 

and prove equations (j27|) . 

For i?^,.,, = ^ - + Li,C;, - - L^,q,J + C}A from 

(|A.6p . we obtain zero values because C**^ = and L*^ do not depend on y'^. 
So, Rja = R'jia = 0- ^ 

Taking i?'^,,, = ^ _ ^ + C1C«,, - C^.C.'^,, from ([MD and^con- 

tracting the indices in order to obtain the Ricci coefficients, R^c = — 

^ + a^^^Ce - d^^di, we compute i?,, = (C^J* + (5,^^° - 5,5, + C^^Ca + 
^^^^^4 - C^gCf^ - C3jC|, - C^^df^ - Cl^C\^. There are nontrivial values, 



4 

34 "-"43 



i?33 — {Cz'i) + (^33) ^ C*! + C'IsCa + CI3C4 — CI3CI3 — 2CI4CI3 — CI4C, 
/i* \ * (hi h\\* hi ( hi h\\ ( hlV f hl^^ 



2/13/ V2/i3 2/14/ 2/i3 V2/i3 2/14/ v2/i3/ V2/i4 
1 h%* 1 1 hi hi 

"2 7^7 ^ 4(/i4)2 ^ 4/^7^1' 

(/-i3 \ I / /-i4 \ fi _L_ I /-i4 /^3 /^3 0/^3 i^4 i^4 /^4 

^44 ) ' 1 ^44 ) ~ C4O4 ~t O44O3 -f- 04404 — O43L/34 — ^044^34 — O44O44 

1 (h\\* h^ hlh^ f hi hi hi' 



2 \hz J h^ 2/13/13 \2/i3 2/14 ft,4 
27^-4+ 4(7^-4+ 4/^/^-4- 



We get the nontrivial v-coefficients of the Ricci d-tensor, 
i.e. the equations (f25]l . 
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B.4 Zero torsion conditions 

We analyze how to solve the equation 

^4 ?4 /a74n dkih^hi) hi III 

T,, = L,, - e4N,) = - w,— - n,— = 0, 

which follows from formulas (jB.Sp for a vanishing torsion for D. Taking any 
/14 for which 

the condition n^-^^ — W^^t^ = is satisfied. For instance, parametriz- 
ing ZL4 = ^lhi{x^)!l{y'^)j the equation ()B.6P is solved for any 

h{y^) = e^y^ and = >cdk[ ^h-i{x^)\^ for x = cont. 

We conclude that for any and h^^ related by conditions (IB.6P the zero 
torsion conditions (jB.SP are the same as for = const. Using a similar 
proof from [ini E] , it is possible to verify by straightforward computations 
that = if the equations ([15]) are solved. 

B.5 Geometric data for diagonal EYMH configurations 



The diagonal ansatz for generating solutions of the system (j30p - (j32p can 
be written in the form 

°g = ° gi{x^)dx' dx' + °ha{x\x^)dy'' (^dy" = (B.7) 
= q~^{r)dr O dr + r'^d9 ® dO + r'^ s\v? 6dip dip - a'^ {r)q{r)dt ® dt, 

where the coordinates and metric coefficients are parametrized, respectively, 

= (^a:' =r,x^ = e,y^ = ^,y^ =t), 
°gi = q-\r), °g2 = r\ °hs = sin^ 9, °h, = -a\r)q{r) 

for q{r) = 1— 2m{r)/r — Ar^/3, where A is a cosmological constant. The 
function m{r) is usually interpreted as the total mass-energy within the 
radius r which for m{r) = defines an empty de Sitter, dS, space written in 
a static coordinate system with a cosmological horizon at r = Tc = y^S/A. 
The solution of Yang-Mills equations ([30]) associated to the quadratic metric 
element ()B.7p is defined by a single magnetic potential a;(r), 

°A = °A2dx^ + °A3(iy3 = — [uj{r)Tide + (cos^ rg + w(r)r2 sin0) d^] , 

2e 

(B.8) 
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where ti,T2,T3 are Pauli matrices. The corresponding solution of ([32]) is 
given by 

$ = °$ = ro(r)r3. (B.9) 

ExpUcit values for the functions a{r),q{r),uj{r),w{r) have been found, for 
instance, in Ref. [18] following certain considerations that the data (|B.7p . 
(jB.SP and ()B.9p . i.e. [ °g(r), °A{r), °^{r)] , define physical solutions with 
diagonal metrics depending only on radial coordinate. A well known diag- 
onal Schwarzschild-de Sitter solution of ([30 |) - ([32]) is that given by data 

w(r) = ±1, o-(r) = 1, (p{r) = 0, q{r) = 1 - 2M/r - Ar^/3 

which defines a black hole configuration inside a cosmological horizon be- 
cause q{r) = has two positive solutions and M < 1/3\/A. 
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